


Institutional Archive of the Naval Postgraduate School 


Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


1986-12 


Remote sensing of ocean sediment volume reverberation 


Chang, Chin-wen 


http://ndl.handle.net/10945/22072 


Downloaded from NPS Archive: Calhoun 


| Calhoun is the Naval Postgraduate School's public access digital repository for 
D U DLEY research materials and institutional publications created by the NPS community. 
3 Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 
KN OX appointed — and published — scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 
Monterey, California USA 93943 





http://www.nps.edu/library 


DUDLEY KIO LS SOARY 
NAVAL bey - SCHOOL 
MON fEREY, Unsikunwla 93943-5002 

















NAVAL POSTGRADUATE SCHOOL 


Monterey, Galifornia 





THESIS 


REVO TESoehoiiG Or OCEAN SEDIMENT 
VOLUME REVERBERATION 


by 


Chang, Chin-Wen 


December 1986 


Co-Advisor James V. Sanders 
Co-Advisor Alan B. Coppens 





Approved for public release; distribution is unlimited. 


7230165 





- _—_ Sa 





SECURITY CLASSIFICATION OF THiS PAGE 
REPORT 


13 REPORT SECURITY CLASSIFICATION 


23 SECURITY CLASSIFICATION AUTHORITY 
UNELASSTEIED 
2d DECLASSIFICATION / DOWNGRADING SCHEDULE 


4 PERFORMING ORGANIZATION REPORT NUMBER(S) 
ba NAME OF PERFORMING ORGANIZATION 


Naval Postgraduate School 
65 


& ANDRESS (City State and ZIP Code) 
Monterey, California 93943-5000 


B83 NAME OF FUNDING/SPONSORING 
ORGANIZATION 


Bc ADDRESS (City State and ZIP Code) 


11 TITLE (include Security Classification) 


Eb OFFICE SYABOL da 
(if applicable) 


Bb OFFICE SYMBOL 
(if applicable) 


DOCUMENTATION PAGE 


Ib RESTRICTIVE MARKINGS 


J DISTRIBUTION AVAILABILITY OF REPORT 
Approved for public release; 
distribution unlimited 


S$ MONITORING ORGANIZATION REPORT NUMBER(S) 


N&ME OF MONITORING ORGANIZATION 
Naval Postgraduate School 


To ADDRESS (City State. and ZIP Code) 


Monterey, California 93943-5000 







9 PROCUREMENT INSTRUMENT IDENTIFICATION MtiMaAER 


10 SOURCE OF FUNDING NUMBERS 


PROGRAM PROJECT 
ELEMENT NO NO) 


TASK 
NO 


WORK whit 
ACCESSION: fet) 





REMOTE SENSING OF OCEAN SEDIMENT VOLUME REVERBERATION 





+2) PERSONAL ALITHOR(S) 
Chang, Chin-Wen 

“Ja TyPe OF REPORT 
Master's Thesis 

"6 SLPPLEMENTARY MOTATION 


13m TiAE COVERED 


FROM TO 


ey COSATi COGES 
FELD GROUP 


eee 





——o—— 


SU8 GROUP 


——— 


IS) PAGE COUR 


118 


14 DATE OF REPORT (Year Month Day) 


986 ebe 


18 SUBJECT TERMS (Continue on reverse if necessary and identify by block number) 
volume reverberation, waveform envelope, transducer, 
receiver, sound speed, slope, decay constant, 


Scqttecina cOehr Ie Lent 


"9 ABSTRACT (Continue on reverse if necessary and identify by block number) 


An experimental study was performed to establish a technique, for neasueing phe voluty 
Bae ey tt ion from ocean sedinents. Two types of sediments (aggregate and fine san 5 


were used in this study. The inhomogeneity within the sediment caused considerable 
sample-to sample fluctuation in the scattered waveform. This fluctuation was removed 
by spatial averaging to obtain a mean value over the sampling area. A approximate 
nodel for volume reverberation from sediments was developed for an acoustic pulse 
with an exponential decay. The results are promising. Combining the model and 
experimental results, the volume scattering coefficient obtained for the aggregate 
is 0.0624 40.007 m and for fine sand is 0.0413 0.007 m The latter result 1s 


NNN 
£0 OS R:BUTION/ AVAILABILITY OF ABSTRACT 21 ABSTRACT SECURITY CLASSIFICATION 
Cy) CLASSIFIEOUNLIMITED C} SAME AS RPT () OTIC USERS ° ied 
229 NAME OF RESPONSIBLE INDIVIDUAL 22b TELEPHONE (include Area Code) | 22c OFFICE SYMBOL 
James V. Sanders/Alan B. Coppens 408) 646-2117 \Sd 


BIAPRedition may be used untlexhsusted 


DD FORM 1473, 84mar SECURITY CLASSIFICATION OF THIS PACE 


Allother editions at obsolete 


SECURITY CLASSIFICATION OF THIS PAGE (Yhen Data Zntered) 





19. ABSTRACT (continued) 


close to the coefficient cbtained from a cloud-or-cma eee 
spheres model for wavelengths much greater than the 
particle size. 


Cee ORE eT, as 
A Ty eT a en ED GELATO cl a SOTA A ET Os cis. A ire es Se ARN fe A st i UR A eee apt Estes 


S/N 0102- LF- 014-6601] 
Z 


LS rr hf phen 
SECURITY CLASSIFICATION OF THIS PAGE(Whon Data Entered) 


Approved for public release; distribution is unlinited. 


Remote Sensing of Ocean Sediment 
Volume Reverberation 


DY 


Chang, Chin- Wen 
Lieutenant, Republic of China Navy 
B.S., Chung Cheng Institute of Technology, 1979 


Submitted tn partial fulfillment of the 
requirements for the degree of 


Mee hel SCIENCE INS ROGRAPTIIC SCIENCE 


from the 


NAVAL TOSTGRADUARE SCHOOL 
December 1986 


Ps oi) penne 


An experimental study was performed to establish a technique for measuring the 
volume reverberation from ocean sediments. Two types of sediments (aggregate and 
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result 1s close to the coefficient obtained from a cloud-of-small-spheres model for 
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lL INTRODUCTION 


In hydrography, the character of the bottom should be determined for nautical 
charting, particulariy in harbors, anchorage areas or other areas where the bottom 
characteristic is a significant factor to be considered for navigation safety. 
Conventional sampling methods, using clamshell snappers (Umbach, 1976). are used to 
define the characteristics of the bottom surface laver. It is tedious and time 
consuming, especially in off-shore geologic investigations. If some kind of theory or 
Miodeme sisted to remotely classify the sediments, the task could be carried out easier 
and faster. Moreover, this ability would permit nautical charts to contain much more 
information. 

Poerie Oceam OOiom, the magnitude and mature of the sound scattering are 
functions of both the particle size of the dottom and surface (bottom) relief (\{ckinney 
peeeicersor, $96.5 Tie sound echo ‘rom sediments ts subject to three tvpes of 
felecuiie Or scattering (Clark, Pront, Seem, Tsai. 1984): 

1) coherent reflection from the water-sediment interface; the signal amplitude 
will be approximately that given by spherical spreading of the wave front 
from the image source modified by the plane wave reflection coefficient 
(G@iay, 1966); 

2) incoherent or statistically variable scattering from the surface irregularities 
of the sed:ment, and 

3) incoherent scattering from within the volume of the sediment caused by the 
acoustically irregular matrix of the sediment. 

Items 1 and 2 arise from the surface of the sediments and item 3 from the sound 
pulse penetrating into the sediment volume and then being scattered back to the 
transducer by the granular matrix of the sediment. Since different sediments have 
different acoustic matrices, they have different volume reverberation characteristics and 
there should exist a mathematical relationship between sediment properties and the 
volume reverberation. 

The major goal of this experiment was to establish a fundamental laboratory 
procedure for the remote sensing of ocean sediment volume reverberation. This was to 


be accomplished by accurately measuring the properties of a pulse reflected from the 
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sediment. In application, since the reflected sound pulse is rece eGupe tices a. 
transducer, Le., “monostatic sonar’, our discussion ts restricted to “back-reflection” and 


“backscattering”. 
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Hl. BACKGROUND 


Bradshaw (1981) used Monterev #31 fine sand, the same material used in this 
experiment, to study the sound propagation into a fast bottom medium. In his study, 
air bubbles were removed from within the sediment based on heating the water sand 
system. Bleach was added at a ratio of one-half gallon bleach to 70 gallons of water 
plus sand in order to control biologic growth. The physical properties, such as density, 
porosity, speed of sound, reflecton coefficient etc., Were measured. Kosnik (1984) used 
a high pressure water jet to remove the air bubbles from the water) sand inixture. 
However, no volume scattering measurements were made by either investigator. 

Morse and Ingard (1968) give a theoretical derivation on incoherent scattering. 
Under the assumptions that : 1) the scattering objects are quite smail compared with 
wave length A. 2) the dimensions of scattering region are niuch larger than A, 3) the 
scatterers are all spheres, and 4) the scatierers are populated sparsely and no inuuple 
scattering existed, the incoherently scattered intensity 1... per unit meident intensity 1, 


Smale g to the incident wave, for frequency 0/27, is 


Lol =NV any? KA a8 fy. +p cos6|? (2.4) 


LK? a2 sin?(1/2)0) exp(-2K? a~ sin?(1/2)9) 
8 


where N is the total nurnber of scatterers within unit volume, K is the wave number for 
the mixture, a is the radius of scatterers. Furtherniore, for monostatic sonar, the 


volume scattering coefficient for unit volume can be reduced as 


sy = N (2x)! ?KS aS] y,. - yf? (2.2) 
munere 
Kn = Ko 
Ye 
o 


I 
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and K, 1s the compressibuity of the scatterer, K, is the compressibility of the fluid. p,, 
is the density of the scatterer and 9) is the density of fluid. 
For a sinall (ka << 1) rigid Sphere, the Sackscatterine ciosc seciiG mise 
(UniCkKe sats 
9 Ss 
G = 2.8(ma“\(K aj" 
Note that, as a single sphere. it must be a single scattering. The scattering coefficient 


for unit volume can be obtained as 


Both theories are not very appropriate for the rea! ocean sediment voluine 
reverberation. Llowever, these are the only two theonmes avatlables Une backseatvenm: 
of sound from the ocean bottom has been studied for many years, but unfortunately no 


papers were found to discuss the volume reverberation of ocean sediments. 
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Nil. DESCRIPTION OF VOLUME REVERBERATION 


As sound propagates in a fluid. part of the sound is intercepted and reradiated by 
inhomogeneites in the uid. The reradiation of sound is called scattering. The total of 
memec elias COMmOUtONS from all the scatterers is called reverberation (Urick. 1955). 

\ olume reverberation is defined as sound scattered back to the transducer bv 
Scattering centers in the volume of the sea (National Defense Research Commitee, 
1969). Volume reverberation mav be considered as a blending of large number of 
pemocseeit Mas been studied mostly in the Deep Scattering laver , a region of laugh 
bieroeical activity. Reverberation level, in decibels, is used to measure yolume 
reverberation. Consider a directional projector in an ocean containing a large number 
of volume scatterers (Figure 3.1). Let the axial intensity at unit distance be 1,. Ata 
distance r where there is a volume element dV of volume scatterers, the incident 


Pena 1s 


Mie scattered intensity { m from the volume dV 1s 


di, (1) = sy d¥ 


“4 


where s. is the volume scattering coefficient. The scattered intensity back at the 


= ae 
eee) = Ls. - 
el Aus ay. 


and the scattered intensity from all scatterers contributing to the scattered intensity at 


time t is 
I(r) = r?) Ss. V 
S O vi 


where V is the total volume contributing for the geometry of Figure 3.1 calculated by 


dieiex, 1953) 


ao r? Cr. 
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TRANSDUCER 





Figure 3.1 Portion of ocean scattering. 


IS 


where & is the solid angle of the sound beam, T is the pulse length and c is the sound 


Speed, so 


) 


Pee) — Cibmi Sena cc, 2 


0 


Tne reverberation level is 


RL = 10 LOG I,{r):Lor 
= SL-TL +S, + 10 LOG (8 ct-2) 


Where TL is the transmission loss for spherical spreading 20LOG(r) and Se ile 
scaitering strength for unit volume 1OLOG(s,.). In Deep Scattering Layer studies, this 
reverberation level equation ts sufficient, because the received pulse is regarded as 
totally conung from: volume reverberation, i.e. negligible reflection and surface 
scattering. Furthermore, the derivation of V is based on a square acoustic pulse. But 
for a real transducer, the square pulse is transformed into a pulse with an exponential 
rise, a flattened top and an exponential decay. The sound pulse reflected from the 
sediment and received by the transducer combines coherent reflected. incoherent 
scattering from the surface and volume reverberation. | 

Let a transducer located in the water send out a sound pulse to a smooth 
Petemmiiitciiace. [he received eclio waviciormi 1s purely reflection. If the water 
interface 1s rough, then the received echo waveform includes both reflection and surface 
Seattemme. ince SOund energy is scattered by the rough surface, the received waveform 
amplitude 1s reduced. Consider first-order scattering and suppose that the surface 
roughness is the same for every position. The ainplitude of the combined sound wave 
(reflected and scattered) 1s reduced proportional to the purely reflected wave amplitude. 
In Figure 3.2(a), the solid line represents an ensemble of reflected waveforin, and the 
dashed line represents an ensemble of reflected and scattered waveform. 

The contributions to a pulse reflected froin a sediment surface is shown on 
Figure 3.2(6). The dashed line is the contribution by reflection and surface scattering: 
again. only first order scattering 1s considered. The dotted line includes volume 
reverberation. The level of the signal is raised above that of a dashed line because of 
energy returned from within the sediment. In particular, the tai! decays more slowly 
than for reflection from the water/air interface. This is the significant feature of 


sediment volume reverberation that we are after in this paper. 


io 





WAVEFORM ENVELOP FOR WATER/AIR SURFACE 





WAVEFORM ENVELOP FOR SEDIMENT SURFACE 


Figure 3.2. Theoretical echo waveform for water and sediment. 
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Two alternative assumptions have to be made for reflection from the sediment 
surface : either !) the surface scattering is negligible compared with volume 
reverberation, or 2) if it is not negligible, it must have the same amplitude factor as 
memsced im treure 3.2(b). 

According to these assumptions, if a sediment echo is subtracted from the water 
Eonemtr il Sive a measure of the volume reverberation in the sediment (Ficure 3.3). 
mime cecaviiie tail of the echo reflected irom: the smooth water, air interface is an 


exponential decay. 


where V,..(t) is the received voltage at time t. 6 is the decay constant. Vi, 1s a constant 
voltage, and t 1s the time starting from the beginning of decay. If we assume the 


sediment echo also decavs exponentially, then we will obtain a simular function 
ee we. ct 

Where V.(t) is the received voltage from sediment at time t and 4 is a decay constant. 

If V,,.and V, are both normalized to the same initial value k, the difference becomes 


CU sas fenle= et) (ely) 


If both yt and ot are much less than unity, 
eae aa = Oe fac Paand Io 


This initial decay is linear and depends on the difference 6-1. 
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A. MATERIAL AND TANK SELECTED 

Fresh (tap) water, #30 Monterev fine sand and Monterey Aquarium #2 fora 
rough gravel sediment were the media used in the experiment. The grain size of #30 
fine sand varies from 0.15 mm to 0.7 mm (Bradshaw, 1981) with an average size 0.3 
mm. The Aquarium #2 has a mean size 5.3 mm (Diaz. 1986). The sand' water mixture 
Was repeatedly stirred to remove air bubble trapped while transferring the sand to_ 
different tanks. 

Initial experiments were conducted in the NPS anechoic tank to measure the 
reflection from the water/air surface. with the transducer located at different depth. 
The second tank, a wooden tank measuring 80 cm x 80 cm x 60 cm, was constructed to 
hold about 50 cm of aggregate. This tank was lowered and suspended 85 cm below the 
water air surface of the anechoic tank. A third tank, a steel glass tank, measuring 70 
cm x 7/0 cm x 60 cm was filled with 24 em of fine sand or aggregate with 36 cin of 


water above. 


Poe bec TRONIC EQUIPMENT 

A schematic drawing of the equipment is shown in Figure 4.1. Except a high 
pass filter, all components were off-the-shelf. General Radio model 1310 oscillator 
with a frequency of 182.0 kHz was fed simultaneously into a Hewlett-Packard 3233L 
frequency counter, and a Gencral Radio Type 396-A tone burst generator was used to 
generate a 16 cycle pulse. Signal output from Ilewlett-Packard IIP 467-A amplifier 
was fed through a Datasonic Transmit, Receive (T,;R) switch and a small 2] to 28 kIlz 
high pass filter, then into the F-4] transducer. The F-41 transducer was discussed in 
detail by both Borchardt (1985) and Diaz (1986). Signals returned from the sediment 
Were amplified 20 dB by a Hewlett-Packard 465-A pre-amplifier, then passed through a 
Spencer-Kennedy Laboratories, Inc. model 302 variable electronic filter (set at 135 KI1z 
high pass) to eliminate low frequency noise before being passed to the digital and 
analog oscilloscopes. 

The waveform displayed on 3091 Nicolet Digital Oscilloscope was sampled and 
sent as a stream of characters to the computer. To reduce manual operations on the 


Nicolet 309! during the experiment, a IIP-3421A Data Acquisition and Control Unit 
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Figure 4.1 Schematic of equipment. 
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was adopted to recerve commands from a HP-86 computer and then automatically 
Operate the Nicolet 091. However, it was noticed that if the Nicolet 3091 was not 
warmed up” enough, the DC-offset jumped up and down, and continua! manual 


adjustment was needed. After 2 or 3 hours of running. it became stable and no further 


adjustinent was require 


See ERIFICATION OF MEASUREMENTS 

No matter how deep the transducer is located below the water air surface, as 
long as the transducer is level, the exponential decay of the tail should be the same. In 
Diaz's thesis, the measurements Were inconclusive because the transducer was fixed at 
one position throughout the experiment. This was not adequate to prove the system 
functioned properly. To prove that the previous experiments were run properly and 
results were obtained correctly, a verification was attempted by measuring the 
reflection from the water, air interface with the same procedure, but with the transducer 
eee Ovo tis, Iain set 20 (Figure 4.2) for the transducer 70 cm below the suriace, 
ees a S.Cpe Of -3.°3 x 07 Np us. Data set 2! (Pisure 4.3) was obtained ov raising 
tie iransdticer 15 cm and the slope is -3.77 x 10" Niece GUM souls Were ciagserG 
MmeeelteviOls Tesuits obtained bv Diaz (1986). After these two data sets were collected, 
a new transducer niounting was built and adopted to improve the experimental 
precision and operation ability. The transducer was leveled by a counter-werght on the 
new carriage. To ensure the system) and operation worked properly, more locations 
[cerned (oO Obtain the decay constant of the water, air surface. 

Data sets 23 and 24 were obtained at the same transducer depth, but seperated 
40 cm horizontally. Data sets 24 and 25 were obtained at the same horizontal position 
but with 50 cm difference in depth. Data set 26 was acquired by changmg 10 cm mn 
Seumtne dept) and position from data set 25. Data sets 26 to 28 resulted from 
eieneame the depth from 70 cm to 30 cn below the water/air interface. The results of 
these experiments are listed in Table | and show consistencies between each other. 
Waveform enseinbles of the tail for each data set are plotted on Figure 4.4 to Figure 
4.9, 


7 - -2~- ; 2 
The average slope -3.77 + 0.077 x 10°* Np. {ts 1s smaller than -3.84 = 0.14 x 
10°“ Np/ pts which was measured in Diaz’s thesis. Statistically, these two slopes are the 
same. The smaller standard deviation nught suggest that the new carriage has achieved 
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Figure 4.2) Decay envelope for reflection from water/air interface, 
transducer 70 cin below the surlace (data sect 20) 
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Figure 4.3 Decay envelope for reflection from water/air interface, 
frinscdneer sscti oclow Wie cupiaceddata set 21). 


TABLE 1 
WATERVAIR Rie Peerie.~. 


DATA SLOPE CORRELATION 
SET (x1074) COEFFICIENT 
DATA20 aes -0.9984 
DATA21 =a “0. O7c5 
DATA22 Some -0.9979 
DATA23 ee ee 2009774 
DATA24 -3. 68 -0. 9895 
DATA25 -3. 86 -0. 9956 
DATA26 =auco -0. 9679 
DATA27 -3.84 -0. 9888 
DATA28 aan 7o -0. 9810 
MEAN SLOPE = ~3.77 x 1077 Np/is 


¢ = 0.07 x 1074 Np/pts 


After dealing with the water/air surface, the transducer was turned over and 
pointed downward toward the sediment which was contained within the 80 cm x 80 cm 
x 60 cm wooden tray. Different transducer positions were choosen. Data sets 31 to 33 
were at the same depth, but horizontally 3 cm awav from each other. The results show 
the waveforin were totally inconsistant. Data sets 33 to 35 were at the same horizontal 
position but different depths; results were also inconsistant. Above results are 
presented on Table 2 and plotted on Figure 4.10 through Figure 4.14. So far, 
experiments reveal the lack of isotropic homogeniety of the volume contributes 


reverberation. 


D. MEASUREMENT CONSIDERATION 
During the course of these experiments some other effects were considered and 


studied to ensure that the experiment was working properly. 
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Figure 4.4 Decav envelope for reflection from water/air interface, 
transducer was 30 cm below surface (data set 23). 
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Pigure 4.5 ] 
trnsducer was 30 cm below surface (data set 24). 
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Decay envelope for reflection from water/air interface, 
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Figure 4.6 Decay envelope for reflection from water/air interface, 
transducer was 80 cm below surface (data set 25). 
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igure 4.7 Decay ens for reflection from water/air interface, 
transducer was 70 cm below surface (data set 26). 
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Figure 4.8 Decav wea for reflection from water/air interface, 
transducer was 50 cm below surface (data set 27). 
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Pigure 4.9 Decay Cree for reflection from water/air interface, 
transducer was 30 cm below surface (data set 28). 


TABLE 
Sew. beaver Weerree MON sm sCATTERING 


DATA SLOPE CORRELATION 
SET Gator) COEFFICIENT 
DATA31 ee Seyoise 
BETAS 2 aye -0. 8908 
DATA33 =o -0. 7838 
DATA34 -3.06 Oued 
DATA35 2590) - >: -0.9844 


MEAN SLOPE = -2.62 x 107% Np/ts 


¢ = 0.86 x 1074 Np/ns 


1. Effects of Sampling Rate 

The sound signal being used was fixed at 182 kIIz (5.5 pts period). The 
Nicolet 3091 digital oscilloscope has as its fastest sampling rate | ts which is 
inadequate to give a precise reproduction of the numeric waveform. However, 
observations have shown “jitter” in the triggering of the scope, so for consecutive 
pulses the scope would sample different parts of the waveform. To obtain an accurate 
waveform, a large number of waveform were sampled and averaged. If sufficient 
Sievers 1s present, tlre voltage measured in a stven btn should vary between the 
maximum and minimum voltage in the waveform in the vicinity of the corresponding 
time. Since the sine wave varies between minus and plus, the voltage values were : 
squared before averaging in order to give an accurate representation of the waveform 
envelope. 

To determine the maximum number of samples required for each data set 
(average waveforin), the water’air interface was measured with different number of 
samples. The results showed that 50 samples for each data set were enough to give a 
true waveform. 

2. Effects of Bottom Reflection Interference 
For a 16 cycle signal with a frequency 182 kIIz, the pulse length 1s 88 Hs. 


Since the sound pulse emitted had a decay time of approximately 100 pts, the echo from 
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Decay envelope for reflection from aggregate/water 
interface, transducer was 30 cm above surface (data set 31). 
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Figure 4.12) Decay envelope for reflection from aggregate; water 
Wnteiface, transducer was 30 cm above surface (data set 33). 
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Figure 4.13 Decay envelope for reflection from aggregate/water 
interlace, transducer was 50 cm above surface (data set 34). 
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Figure 4.14 Decay envelope for reflection from aggregate/water 


interface, transducer was 70 cm above surface (data set 35). 
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the smooth water air mterface should also have the same decay time. Since the 
purpose of this project was to compare the decay portion of echos reflected from water 
and sediments, the waveform was sanipled for a 100 js interval from the begining of 
decay. The sound speed tn the aggregate was measured to be 1555 m's. For the fine 
sand. Bradshaw (1981) obtained the speed of sound to be 1619 ms. Combining these 
sound speeds and the time interval 188 us, the depth of sediment that contrubutes to 
iewmiat cco 1$ 1456 ci for acsrecate and 13.1 em fer fine sand. Interference caused 
bv refiection from the bottom of the tank will not affect the experiment as long as the 
sand is thicker than 15.1 cm and the aggregate ts thicker than 14.6 cm. Throughout 
Our experiments, We kept both sediments more than 24 cm deep tn the glass tank to 
avoid mterference. 
3. Effects of Near Field Interference 

Mimiicenedy eldeo! 4 transiiutter, te axial pressure eXhibits strong interference 

Biteerse hetisier, Frey, Coppens., Sanders, {9S2). fhus. the sound field is irregular and 


does not fall off smocthly with distance as it does rn the far held. To avoid 


iem@eseoarates the near field and far field. This distance ts 20 cm from the transducer 
center. No problems should occur as long as the target ts more than 20 cm from the 
transducer. In the glass tank, we placed the transducer 30 cm above sediment. The 
Water air interface measurement also gave a check that 30 cm distance 1s satisfred. 
Since a pulse moving in water at 1480 m’s, takes 405 [ts to travel the 30 cm to the 
sediment interface and back. the 30 cm distance ts long enough so that the transducer, 
which is on about ISS Its, stops ringing before the first echo from the sediment arrives. 
+. Effects of Water/Air Interface Interference 

Ideallv, the transducer should be far enough away from all unwanted surfaces 
to avoid interference. As the transducer faces down and projects a major part of Its 
energy to the sediment surface, some energy 1S transmutted to the water surface and 
reflected back toward the sediment. If the transducer ts close to the water, air surface, 
the reflected energy will interfere with the transmitted energy. Restricted by the size of 
the glass tank, we had to locate the transducer close to the water surface. To test how 
Iniportant this interference could be, the water surface was agitated by fingers. The 
waveform displayed on both analog and digital oscilloscope was unchanged during the 
ruffling. It was concluded that this source of interference is neghtgrble for thrs 


exe wMent. 


4} 


5. Effects of Side Lobe 

The half beamwidth (to the first null) for the F-4} transducer was measured to 
be 10.0°. For a 30 cm distance from sediment to the transducer, the sediment suriac 
ensonifted bv the major lobe is a circle area with radius 5.3 cm. If the distance is SO 
cn. the radius will be 14.1 cm. For both the glass tank and wood tank, if the 
transducer is located above the center of the tank, tie major lobe will moc srocuce 
interference caused Ob. tie sieec te) We ts 

The glass tank was too small to avoid the interference, if anv, caused bv the 
side lobe. To test how much the interference was, secliments outside (ie cenicn- Ure miee 
30 cn) square area Were disturbed and left in a very irregular surface. Waveforms 
returned from sediment were measured before and after the disturbance. Data set 51 
(Figure 4.15) shows the waveform measured before the sediment was disturbed and 
data set 52 (Figure 4.16) shows the waveform after disturbance. Bv comparing the two 
waveforms, it is found that the F-41 trandsucer gave negligible side lobe interference in 
this experiment. 

6. Effects of Sediment Inhomogeneity 

Before the test findings on side lobe effects, we were bothered bv the drastic 
changes on the waveform whenever the transducer was moved parallel to the smoothed 
aggregate surface. It was assumed that the waveform fluctuation was caused by the 
inhomogeneities of the aggregate. To prove this, we placed the transducer at a fixed 
position in the glass tank and repeatedly ensonified the aggregate. After each data set 
was acquired, the aggregate was disturbed and smoothed again. Ten sets of waveform 
data Were obtamed and each of them differed from the others. 

It was thought the fine sand would be more homogeneous than aggregate and 
it mught produce more consistant waveforms as the transducer was moved horizontally. 
\leasurements made with the fine sand, however, showed that the waveform still 
fluctuated. After each data set was obtained from the smoothed fine sand, the 
transducer was moved horizontally to different position. All the transducer positions 
were at least 15 cm from the wall to avoid interference caused by reflection from the 
sides of tank. Also, the positions were at least a length of tranducer diameter from 
each other to avoid a waveform replication. A metal plate was set on top of the sand 
and the transducer was moved horizontally. As expected, the waveform on the 


oscilloscope was unchanged but totally different from sand. 
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Figure 4.15 Decay envelope for reflection from a 
interface in the glass tank (data det 5 
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Figure 4.16 Decay envelope for reflection from aggregate, water 
interface in the glass tank (data set 


E. PROGRAMMING 

‘An HP-S6 desktop computer was adopted to collect and analyze the waveform 
data. The computer program) was written to be interactive. This enables 1) the 
collection cf data sets with different delay time intervals and 2) the collecting and 
ns of numerous data sets. The interactive programs allowed essential 
information to be entered via the key board. eliminating the need to modify the 
program whenever data sets were changed. The interactive prompts were useful for 
people unfamiliar with Basic programming. 

eeroetame, Cie SIS 5” 

This program was designed as the initial step in the data collection processing. 
Before the selected waveform data displayed on the 3091 Nicolet Digital Oscilloscope 
MenesOuUtpul TO [he COMpluter, ii was switched into “freeze” or ‘store’ mode from live” 
eee) nroush the mterface RS-232, the conipttter received the stored waveform data 
and then processed them. After processing the waveform data and before receiving the 
iecwsaiapume data, the \tcolet 35091 had to be switched into ‘hve’ mode again. An 
HP-3-21A Data acquisition and Control Unit was interfaced between the Nicolet 309! 
and HP-86. The HP-3421A reduced the amount of manual operations in the data 

ollection process. A five volts source were input to the HP-3421A. Commands form 
the HP-86 and 5 volts signal were used to control, by closing and opening the channel 
in the HP-3421A, the “mode switching” of the oscilloscope. To avoid errors, the 
oscilloscope should be warmed up enough while entering the imital information. Once 
this initial information was entered and the oscilloscope was warmed up enough: the 
Whole svstem runs automatically. 

Waveform data output from the oscilloscope were a stream of characters. [To 
obtain the voltage and time vatue, all these characters meeded to be normahzed and 
converted by the computer. A program, listed in the Nicolet 3091 operation manual, 
was used for data value normalization. A data set, the waveform of the decay portion, 
was obtained by squaring and averaging 50 samples of waveform data for each 


transducer position. The running average formula used in the program was 


ee a en 


2 
Where < V~ >, isn samples averaged squared- voltage, < v2 > n-| 1S n-1 samples 
averaged squared-voltage and ee is squared voltage for the nth sample. Note that. 
the average window for the waveform is 100 [ts from the beginning of the decay. The 


above formula was processed for each corresponding time bin. 
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During each data set processing, some error messages, such as STRING OVE 
(which stands for the string variable overilow)] were displaved on tiie F=sessence 
This error messaged could be fatal to the data set 1f the data has not been stored on 
the disk. To avoid this problem, the program) was designed so that after every 10 
samples processed the updated data were stored on the disk over the previous updated 
data. The >rogram eliminated the problem of the disk being full before the final data 
set was stored on the disk. Tnrts prevention could increase the processing time, but 
compared to the time wasted from the above problems, it is Worthwhile. The final 
result for this progran) was 3 data files : one stored the squared voltage, one for time 
data and one for the number count of the samples. 

2, Program “MIANINIUNI” 

In order to obtained the envelope of the waveform. this program was designed 
to pick up the local maximum? of the voltage data. [lic voltages’ stored aca 
from “TIHESIS3”, were read and compared with the neigh borne voltage datancge eer 
the local maximum. The determined local maxmmiums and corresponding time bins 
Were stored in seperated files anu printed out. [ite cata essere ci crenn encom 
number according to the order they were acquired. 

3. Procritie tee Os 

After the local maximums Were determined, the decay portion of the waveform 
can be plotted by program “PLOT3°. The expemmients, [on thus project, 1001 1ee 
mainly in the glass tank. A decay envelope measured in the glass tank has an interval 
between 520 and 619 [ts starting from transmitting time. Some exceptions existed and 
gave different trme frames. 

To determine the decay constant, the voltage data were converted 


logarithnucally and plotted vs. the time 


LN(V) = -at + LN(V,) 


If the cursor on the oscilloscope was not centered, the resulting envelope showed a 
“Zig-Zag” instead of a straight line for a perfect exponential decay. 


The least-square method was adopted to calculate the best fit of straight line 


y= ax tb 
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Where a, which is similar to the decay constant, is the slope and b is the intercept on 
the v-axis. According to Beers (1953), 


Sree: 
KA (Xp¥p) > 2 Xn 2 Yn 








Nie eNT as . 
Bee a Spb n’ 
—_ 
6 
Oe) 


eeterr ts tie total number Of pairs of (x,y) values, is thle number between | and k, 
a De 
» X,, 1s the sum ot the total x value, ny Y, is the sum of the total y value, v oe es ie 


; oy. 
sum of the squared x value, - Vy 1s the sum of the squared y value and y is 


Lan 


the sunt of the multiples of Xp and ¥p- 


= 


> oe 
prone Chall Ss a and ¥ 2d (X,)> are not the same. y cae MNEs tinat each 


“n 
value of x is squared and then a sum is made of these squares, while ¥ (x,,) lnplies 
that all of the vatues of x added together and then this sum is squared. There is an 

ae Ypy: In our case, the x-coordinate 
is the time value and the y-coordinate is the ack. Oethie wolldce. | Moreneck Now 


analogous distinction between © (Xy¥n) and \) 


close the observed and least-square estimated y-value, compare, a correlation 


coeflicient was calculated (Chatfield, 1984) 


CONG a) 


Oli Gee, se 


where r is the correlation coefficient, x is the incan value of x and vis the mean value 
of y. Slope and correlation coeflicient were both printed with the graph. 
4. Program “AVERAGE” 
The received echo waveform varies from position to position because of the 


inhomogeneous sediment. To find the decay constant, we have to average out the 
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inhomogeneities to obtain a mean value of the waveform. This program performs a 
running average on the original data sets, producing the final averaged waveform of the 
sediment. 


The above four programs are given in Appendix A. 
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Voces tlie 


Mgereata seis for reflection from aesrezate water interface are data 42 to 31. 
frome © -1 to Picure C.9, and Figure 4.15. Table 5 shows the slope obtained for each 


eee 3 setae OE ee Zig = = are =e ay tee : ‘ a toe + 
Mo@desee warieG Cecween Uata sets. [fo Bei a mean value for the sicpe, the running 


so 


Bees emmetaoc was used to process the cata sets. The running average method also 
Sii@ms MOW iast the average processing became stable. Data sess DAT2 to DAT, 
Figure C.10 to Figure C.18, are the results for each average process. The slopes are 
listed on Table 4. It is evident that the slope becomes stable after $ sets of data are 
Seeemaeed.) {he final averged slope is 2.13 x 10"- Np fis and is quite different from the 
Beeuereecced {rom the Water air interiace. The difference between the two slores 
resuited from the elfects of volume reverberation in the aggregate, and the scattering 1s 


Some ait. 


Po lates 
ORIGINAL AGGREGATE REFLECTION 


DATA SLOPE CORRELATION 
SET (x1077) COEFFICIEN 
DATA42 -0.94 Oe 7S 
DATA43 Ste 35 mss 74 
DATA44 -2. 68 -0. 7848 
DATA4S E>. 64 Poueses 
DATA46 Ona BO We? 
DATA47 Sieg encoyin 
DATA48 -0.95 -0. 3990 
DATA9 =o. -0. 8654 
DATASO eos -0. 9756 
DATAS1 0) 28 -0. 0584 


The data sets for reflection from fine-sand water interface are data 55 to 64, 


Figure D.1 to Figure D.10. The slope of each data set is listed on Table 5. Evidently. 
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TABLE 4 
AVERAGED AGGREG. Ie REPEC Tien 


DATA SLOPE CORRELATION 
saz (x1074) CORFEICIENT 
DAT2 Een -0. 6051 
DAT3 ey "Oo 70 
DAT4 Ose =0, Olas 
DATS CS ~0. 9343 
DAT6 iG -0.9914 
DAT7 -1.98 -0. 9952 
DAT a2alc =O, S9ae 
DATS 16 -0. 9943 
DAT10 Sis -0. 9958 


the waveform changes very much for different data sets. The same running average 
method was used to remove the inhoniogentties. Data sets DATI2 to DA] 20 were 
produced from the average process. They are listed on Table 6 and plotted from 
Figure D.1 ito Figure D.19. The slope became stable alter 6 data sets ere, cua 
This might suggest that the fine sand has a greater homogeneity. The averaged slope 
for the fine sand is 2.37 x 1077 Np. pis and is greater than the aggregate. However, the 
Slope difference between them ts not as distinguished when conmipared to the slope 
difference between the sediment and the water? air surface. 

According to Equation 3.1, the decay portion volume scattering, represented in 


voltage, for the aggregate 1s 


V, = K (0213000 ¢737700t y 


and for the fine sand is 
V, = K( a 2oOm 7 e7 3 7700t ) 


The above relations only existed without surface scattering or the surface scattering 


was negligible. 


TABLE S 
ORMGi Nb rie SAND REFLECTION 


DATA SLOPE CORRELATION? 
SET (x1074) COEFFICIENT 
DATASS est aOesige 
DATAS6 Ecos 0, OSenl 
DATAS7 So Ocal? 
DATASS8 -0.75 =O 7290 
DATAS9 OO ome 
DATA6O eA -0. 9957 
DATA61 = 3036 -0. 9050 
DATA62 Se #0. 4920 
DATA63 esl -0. 9389 
DATAS4 ieee -0. 5784 


A simplistic model for sediment volume reverberation was developed and is 
described in Appendix A. Combining the experinental results and the newly developed 


3 for fine sand 


model, the volume scattering coefficient obtained was 0.0413 + 0.007 nv 
and 0.0624 + 0.007 m*? for aggregate. Ihese value are reasonable compared to the 
difference between the decay constant obtained from the experiment. They are 
produced as a best fit using the least square method. The correlation coefficients, 
approximately 0.88, indicate they are quite a good fit. As a first model, the result 1s 
pronusing. Elowever, the discrepancy between the experiment and model nught suggest 
the model is still very crude. Also mentioned in the Appendix A 1s that this 1s onlv a 
first step to try to obtain a theory to describe the volume scattering effect of the 
sediment; it need to be improved bv further effort. 

The ka value for fine sand 1s 0.12, which ts quite small compared to |. From 
Equation 2.2, cloud-of-smiall-spheres model, the volume scattering coefficient for the 


fine sand is 0.0082 m7” and is not too far from 0.0413 n>. 


As mentioned in Chapter 
2, Equation 2.2 was derived under various assumptions and does not apply to real 


sediment. From Equation 2.3, single-small-sphere model, the volume scattering 


2) 


TAB ISE.6 
AVERAGED EINEvs 2 ee ie 


DATA SLOPE CORRELATION 
SET (x1074) COEFFICIENT 
DAT12 eee -0. 9265 
DAT13 aoe -0.9051 
DAT14 =a =O oa05 
DAT15 ecimay =O 79a 
DAT16 23eeG -0.9912 
DAT17 eons -0. 9904 
DAT18 aoEse =0. 99a2 
DAT19 =) se “9, does 
DAT20 “Fea -0. 9916 


“~y 


coefficient is 0.34 m7’ which is an order of magnitude different from the the ineasured 
value 0.0413 m7°. These values suggest that the cloud-of-smiall-spheres model is closer 
to the newly developed model than the single-small-sphere model. 


3 and from 


The scattering coefficient obtained from Equation 2.2 is 9456 nv 
Equation 2.3 is 1134 m™’. Obviously the two values deviate considerably fron) the true 
value. The ka value is equal to 2.1, which 1s outside the Rayleigh scattering limit. 


Thus Equations 2.2 and 2.3 can not be appitediio tlicw teen: 


Vi. CONCLUSION AND RECOMMENDATION 


The following conclusions are possible: 

1. The inhomogenities within the sediment caused the reflected waveform to 
Change considerabiv for various wansdueer positions. The changes can be 
removed by averaging over the inhomogeneities and the averaging yields a 


mean decay constant. 


2, Ine giass tank works very well, and therefore a bigger tank to contain the 
sediment is not necessary. 
3. The laboratory procedure to measure the volume reverberation, using 


off-the-shelf equipment and a high pass filter in a water-over-sediment 
Sesicml. ig WOrka ie. 
4. The basic program, written for a HP-86 desktop computer, is suflicient for 


as eX Deliment. 


Can 


Although Equations 2.2 and 2.3 can not be applied to aggregate since it 
has a ka greater than 1, the aggregate does have a higher volume 
reverberation than fine sand. The aggregate has a larger size than the 
fine sand; this might suggest that the larger size sediment causes more 
slime reverberation. 


6. From the seattering coefficient 0.0413 + 0.007 m7? for the fine sand and 


J 


ee 008m for the agerevate, the scattered energy received by 


transducer is quite small compared to the incident energy. 


—~) 


For fine sand sediment, the cloud-of-small-spheres model works better 
than the single-small-sphere model. 
8. The simplistic model, derived from the true acoustical pulse and described 
in Appendix A, gives a crude formulated volume scattering relationship. 
Combining with the experiment result, the sediment’s volume scattering 
coefficient can be estimated. 
Further study of the volume reverberation from ocean sediment is recommended. 
Future experimentation should inelude further variation of the sediment size and of the 
acoustic frequency to find a relationship combining the ka value and volume 


reverberation. The developed model should be iniproved to obtain a better result 


coincident with the experimental result. The absorption coefficient of different 
sediment with the frequency used should be measured to compare the value obtained 
from Urick’s formuia described in appendix A. Equations 2.2 and 2.3 should be 
improved to obtain the scattering coefficient. Which can be used to check the above 
model. Continued use of the experimental procedure, equipment and computer 


programs, Written for this project on the HP-86, is highly recommended. 
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APPENDIX A 
DERIVATION OF SCATTERING FROM THE SEDIMENT 


The generated electrical pulse can be described by the step unit function. EI, as 


Follows: 
I(t) = ie (et) - tet) ) (A.1) 


where I(t) 1s the intensity at time t, I, 1s the intensity amplitude, t is the pulse length, 


0) 0 
H(t) = 
J tc 0 
and 
0 t-tT < Q 
H(t-t) = 
] t-t > 0 


Since the transducer is a resonant system with daniping, the square electrical pulse was 
transformed into an acoustic pulse with an exponential rise, an approximately flattened 
top and exponential decay. The shape of the acoustic waveform can be described, 


{Sanders and Coppeuis, 1986), as 
for (lee) H(t) + 1 - 8?) et) (A.2) 


where t is the delav time and @ 1s the decay constant. 

Let a transducer, located in the water, send out a sound pulse to a smoothed 
sediment surface. To assure the surface scattering is negligible and the volume 
scattering is only single scattering, reflected pulse received, by the same transducer, 


should have an intensity (Sanders and Coppens, 1986) 


P=.(fT2 se PS 1, mt-28/cydb(r +6)? + R; Mtyir?) / 2? ey 
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Where T; is the intensity transmission coefficient, s,, is the volume scattering 
coefficient, B is the absorption coefficient in the sediment, € is the affected sediment 
depth, R; is the intensity reflection coefficient and r is the distance between the 
transducer and sediment. From above equation, the intensity contributed by scattering 
1S 

2 


25 eo 1, fae creat neta (A.4) 


Pe 


and the intensity: contiiguledspn Wie ve! lection s 


—l 


Ri ty {tr (A.5) 
Ny een 
Te = (ieee ote Wee iit ece wesc (A.6) 
Wile te 

Rt-2&/c) = (1 - oh HE25C) y pye2E cy - (1 - eT AET 25'S) ) H(t 28 ') (A.7) 


During the integration, the lower boundary of € is zero: the upper boundary is ct 2 for 
H(t-25%c) and e(t-t)'2 for H(t-t-2€ ¢).” Note thatonlyhe decay portion.) stam ae 
interested. After the integration, 


1 = A(-(68P'2 - 1) fh GePe 2 co) arene (A.8) 


a (e7CB(t-T) : ae abs (eC B(t-T)/2 b oe ve is TO) 


where 


_ 72 4 
A =e eee (A.9) 


From the experimental formula f}=0.25f (Urick, 1979) and fis 182 kl1z, the absorption 


coefficient 


B = 45.5 dB.m = 5.2 Nom, 


c}.2 is 4030 Np s for aggregate and is 4186 Np’s for the fine sand. The intensity decay 
constant @ 1s twice the pressure decay constant 6. From the experiment, 6 is 37700 
Np.s and it gives an @ equal to 75400 Nps; 2a.c is 97.2 Np’m for the aggregate and is 
93.6 Np m for the fine sand. After piugging in the above values for the aggregate. 


Equation A.$ decomes 


er io) A(813 208) (A.10) 


For the fine sand, Equation A.7 becomes 
ee er. alll) - A(S.6 27") (A.11) 
sy 
BS = R- 
ie = a = | 


Wie woressure reflection cocflicient R is 0.356 for the aggregate (Diaz, 1956) and 0.36 
for the fine sand (Bradshaw, 1981). For the aggregate, 

Reo 1270 

tee 720 


For the fine sand. 


T; = 0.3704 


The puise length is S$ jis and from equation A.5, 


- -OT -75400t 4 nt 
re lamin (e = lene ir (A.12) 


Using the above values, the reflected intensity for the aggregate becomes 
ero OY ir (A.13) 
and for the fine sand 


Weocice OO! | ir" (A.14) 
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The intensity contributed by the scattering for the aggregate is 


1. = (sully fe? 7 970U 1 Gee (A.15) 
and for the fine sand, 
1. = (s, 1, ot" (et I omen cman (A.16) 


Combining the above equations, the reflected intensity from the aggregate received by 


the transducer 1s: 


I = (1, r4)B (A.17) 
where 
B = ((sl@e 0052 Saas 5 ete ore (A.18) 


For the fine sand 


1 = (1,'t7) D (A.19) 
where 
D = (5, ¢ OU Sue Sige og eee a eee (A.20) 
By inputting a series of scattering coefficients and computing the Log (B)! 2 and 


Log (Dy) 2, the theoretical envelope was determined. A least-square-fit slope was 
obtained for each envelope. The scattering coefficient for each sediment can be 
estimated by the value of Sy for which the slope of the theoretical envelope is equal to 
the slope of experiment data. For the fine sand, measured the decay constant 23700 
Np,s gave a scattering coefficient of 0.0413 m°>. For the aggregate, the scattering 


coefficient 0.0624 n> 


was obtained from the measured decay constant 21300 Npvs. 
Figure A. and Figure A.2 show the envelopes determined theoretically from the 


above values of s,, for the aggregate and fine sand respectively. The theoretical 
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cacleme isaplotted Dyenernmalizing the left 


value ee ine exreriment and then 


a 


irom 


cie@wemcal envelore. 


Resulis from the fine sand were sele 


ee 


eo. as Loc 


ad fr 
Nee 


coefiicient cbtaine aceve metnoc. 


os 


Trezli SCeaztcr®e. 


end of the dest fit slope 


~ re = ~ -} - owe 
BE COELICiOn?. Wite Uncemtainsy 


line to the same 


applying this normalization to tne whole 


d to fin O 
By COmparns aysemequol envelopes 
= 9.907 


~ 4s 7 
‘vas G&kctmined. 


- q 4 A «\ ae o as « oe . 2 ae on! Se = “ 2 
and A.3 snow the envelores for scaccering coefficients = 0.007 different from 0.Q-1: 
-3 ! S hances } rer) ° O° sant chanes THe Recast: oe 

me. 1meslope chances as the scattering coefficient changes. In the cegmeing o! 


decay, che dominant term in ane A.29 1s 


~, _ 


on ss and the scattering term whicn 


Tine norme@mzaiion,. de 


ot - esa 
theoretical enveic 


eae c 


° Sareea aN nso rw oe 
eAl oO whe eae 


OR Criss iz 


6] 


mcerimentai enwe! OD : Por ilk 


term which 1s independent 
The theore 


ihe reflection 
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Figure A.| For the aggregate, the decay envelope obtained from the model 
vs. the decay envelope obtained from the experiment. 
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Figure A.2. Vor the fine sand, the decay envelope obtained from the model 
vs. the decav envelope obtained from the experiment. 


63 






LN ¢(<VOLTAGE“2>"1/2) 


SWIL 





= | 
8 120 | 
uh 
7) 
9 
140 Ww in 
2 G&G 
= ia 
a Y 
160 z : 
: & 
180 ‘ 
© 
&3 
200 
220 


Figure A.3 For the fine sand, the decay enyelope obtained from the model 
vs. the decay envelope obtained from the experiment. 
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Figure A.4 For the fine sand, the decay envelope obtained from the model 
vs. the decay envelope obtained from the experment. 
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APPENDIX B 
COMPUTER PROGRAM 


£OOO Bree ee ee eee ee eee eee a a ae ae ee ee 
LOO i 

LOW2ZO Ny 

10S0 Ee fal FROGRAI : THeStS 

1040 thes 

LOS ee. 

160 at THIS FROGRAM IS DESIGNED TO CONVERT AND NORMALIZE THE 
Pose a THE WAVEFORM DATA VALUE AND FOINT LOCATION, TRANSFERRED 
1080 4 FROM THE NICOLLET S091 DIGITAL OSB EEASEDEE-) 10 ene 

1090 a VOLTAGE AND TIME VALUE. THE NIGGEET IS COIROEED TSY 

1100 a4 THE HF-3041 DATA ACQUISITION AND CONTROL UNIT. 

Pio ae THE FROGRAM SQUAFE THE VOLTAGE VALUE AND THEN TAKE 

P20 a THE RUNNING AVERAGE FOR FIFTY SET OF READING. Te 

Lisa U8, INSURE THE DATA IS NOT LOSSED DURING COMFUTER ERROR, 

1140 ae THE DATA IS BEING STORE IN THE DISK AFTER EVERY 10 

Lise? ee SET CALCUEATION. 

160 oe! 

1170 a : 
1180 f PU PERE te ee eee 
f1St 

hie : 


aQOOMSDIM OCLat lo Cl Sl Seles, 

2010 DIM VDATAFEZOI ,VITIMEFl 20a, VCOUNT EC loa 

2020" “INTEGER NGI 

eet : 

2022 % INFUT THE FILE NAME TO STORE DATA 

2) 5 al! 

2<924 DISP “WHAT IS THE FILE NAME FOR VOLTAGE DATA =" 

2025 INFUT VDATA 

2026 DISP “Wheat IS THE FILE NAME EGReiIneE Pata -” 

2027 INPUD RS iihies 

2928 DISF “WHAT IS THe FILE NAME CEeite COUNT Alen. 

2029 INFUT VCOUNTS 

2050 DISF “YOUR INTEREST TIME INTERVAL [1S PROM CSteaTd) remeeibia 
2031 DISF “WHAT IS THE STARTING TIME GF MEASURINGMIN TIERS -SEG> 2] 
2032 INFUT START 

2033 DISF “WHAT IS THE ENDING TIME GE HEASURING (CRIN MICRE-SEE) =~ 
2034 INPUT ENDT 


2083 t 

20.36 ! CHECKING THE INTERESTED TIME INTERVAL FROM CRT SCREEN 

oN ' 

2940 LET N=1 

2060 =! 

2061 ' IF THIS IS NOT THE FIRST SET OF DATA, THE COMPUTER IS GOING 
2062 0.) TO READ THE DISK FOR FREVIOUS DAT INFORMATION, AND CARRY ON 
2065 ! THE, JOB 

2064 { 


2080 DISF "IS THIS YOUR FIRST DATA SE GUT GF THE so Dal eeeHi- 
2100  DISP “PRESS 0) FR Yes oo lee 

2120 INFUT ANS 

2140 IF ANS=1 THEN GOSUB READ _DATA 


2215 
see eles FOLLOWING IS THE COMPUTER RECEIVE THE Dela ERO eo aND 
aoe : THEN CALCULATE THE VOLTAGE AND TIME VALUE. 


39909 i 


22490 DIM I#€5]17,Ds020000] ,NSC40) 
2260 CONTRUMRST eee ors 

2280 CONTROLS os 
22-970 CONTROL ees. 1 
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om -— += mm mee ee —™- = oo fo a 


eae 
Seon 
aad 
ZaOe 
Zous 
2.304 
ESOS 
2306 
2oaus 
2308 
aot) 
eld 
Ne 
21.3 
2314 
Zk) 
2316 
ace 4) 


arn 
fae ne ee 


oss 
nee 


Bae 5 
Zo24 
aoe! 
= 3540 
2260 
ZB 
eee gt) 
D4 
2420 
2441 
ZA41 
2443 
2443 
2446 
PANT 
2440 
=~480 
20} 
Baa) 
one 
235 
25051 
ey eZ 
2540 
2960 
2580 
2609 
246209 
Be2Zt 
2640) 
2669 
2680 
2700 
27 Zu) 
2740 
27460) 
2770 
eee 
2780 
=BOO 


1 
: 
BUTFUT 7O9 s"OFHA4” 


ABOVE STATEMENT IS LET 


FOLLOWING ONE IS TO LET 


t 
: 
{ 
t 
t 
( 
. 
! 
. 
t 


BUTEUT 709 
WAIT 1000 


S~eloo 


Press BOTICH GHLY 


( 

: SED 
Peo Phe SS 

J 

} 

4 


S091, SUTTON” 


LET THE DATA TRASMIT TO THE 


COTE 77 3" Ci-S2" 


J * 
{ ABOVE STATEMENT LET THE WAVEFORM 
1 


Peer G<, 10 COMPUTER THROUGH 


i] 
ENTER 10 3 
ENTER 10 3 
ENTER 10 : 
I 

HO=VAL 
V1= (VAL 
H1=(VAL 
( 
OUTFUT 
OUTPUT 
' LET 
OUTPUT 
aa 
5 
PRINT 


Le 
DF 
tit 


CMrito., 20) 
Wee oie) a) + lOo° (VAL 
fiemes ) 2 oad yo) eit (VAL 


woes OPMS” 
pOoe. OFM” 
THE CART IN LIVE MODE 
POTN CESS 


ie at 


-FOR IT=0 TO ENDT 


MeereownrRt OR IlENDT THEN 2780 


1 
. 


Tee eee 


THE CERT ih 


"NeiIve” “CCHDIT LON 


"STORE" CONDITION 


PITHOUT THe He -s42) UN 


COMFUTER THROUGHR &S-rSe 


INFORMATION BE TRANSFERED 
Ro-7oeo Nlerhreee 


Loe 


(M$026,30])-12) 
(Nf036,403)-12) 


: NORMALIZE THE VOLTAGE VALUE AND TIME FOR THE INTEREST INTERVAL 


D(I)=VAL (DftCJ,KI1)*V1 

T (1) =(1-HO) ¥H1t 

Del -bel) 2 

IF N=1 THEN GOTO 2640 
S(L)=((N-1)*S(1)4+D(1))/N 
I 
IF 
IF 
IF 
IF 


N=1 THEN SCI) =D(I1) 

Pothier tien 2/60 ELSE S680 
Sets 2 HE, 2760 ELSE 

Potent SO) THEN 2760 ELSE 

Peo +SO THEN 2760 ELSE 

Pete ehMbT THEN 2760 ELSE 2780 
eee) Clos SEC Peis VOLT 
mani Wit)? SEC eo Clo Os |. 
' 

J=J+5 

MoV +S 
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2700 
ee) 
Zt A) 


vet By) 
‘a (1) 


aeel) 
2840 
seen 
eo) 
222 
2840) 
ZIe0 
2770 
2971 
Zee 
29780 
27 70 
SOOO 
Sow. 
PRT HY) 
Go7 
Soe 
otis 
Sa 
STOO 
oo i) 
aie) 
2540 
Sog0) 
ao 
S380 
a7 
3600) 
3620 
2640 
36609 
34680 
71) 
reo 
3740 
3740 
7780 
28OD0 
ero 
7 i> 
ey a) 
S750 
3740 
St) 
AOD 
ALO 
AO4ag 
4O560 
4061 
4080 
4100 
4120 
4140 
4160 
4180 
A200) 
4220 
4240 
4260 
4280 
42.80 


MEAT I 
N=N+1 
’ 


: LET THE DATA BEING STORE AFTER EVERY 19 SET OF CALCULATION 


IF N MOD 1O=1 THEM 2945 ELSE 2240 
IF N=11 THEN 2970 ELSE 2980 
CREATE VDATAS, 200,98 
CREATE VCOUNT?,1,8 
CREATE VTIMEF, 200,98 
GOSUB WRITE_DISK 
IF N=Si THEN 3000 ELSE 2240 
END 


: 
' 
: FOLLOWING SUBROUTINE IS USED TO WRITE THE DATA TO 
i 
( 


WRITE _DISK: 

FURGE VYDATAF 

CREATE VDATAS,200,9 
ASSIGNH 1 TO VDATHAE 
FURGE VCOUNT Ft 

CREATE VCOUNTS,1,.8 
ASSIGN¥ 2 TO VCOUNTE 
FURGE VTIMEs 
CREATE VTINEF, 200.8 
ASSIGN# 3 TO VTIMEF 
FRINTH 2: M 

FOR I=START TO ENDT 
PRINTS et (1) 
FRINT# 1 3 S(T) 
NEXT I 

ASSIGN 1! TO »* 
ASSIGNH 2 TO 
ASSIGNH 2 TO -« 
RETURN 


: 
' 
: FOLLOWING SURROUTINE IS USED 10 READ THE DATA FROM 
; WRETTEM Dist 

1 

( 


READ _DATA: 

ASSIGNA 2 TO VCOUNTS 

ASSIGN# 1 TO VDATHAS 

READH 2: N 

FRINT N 

FOR I=START TO ENDT 

READH 1 3 SCI) 

IF I=START THEN 4220 ELSE 4140 

IF I=START+20 THEN 4220 ELSE 4160 
IF I=START+50 THEN 4220 ELSE 4180 
IF I=START+50 THEN 42270 ELSE 4200 
IF I=ENDT THEN 4220 ELSE 4240 
FRINT S(1) 

NEXT I 

ASSIGNH 1 TO * 

ASSIGN# 2 TO »* 


RETURN 
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THESD 1 St. 


PREVIGUS 


Bie ' eres anu enaUe tae mete Noll Pel AY te Onin Me Ae i teak NM Ste, aie ft ak Se AE 
Bl 
Sis 
Sa.3 
B84 | 
SS ! 
BOG 
807 
808 
t 
\ 
t 
t 
' 
j 
' 


Die Lage EC Ee ee inet A Cale tw 8 ne we ia SmE nl Gk be ee ark) costa, Gal” ky Muth Se cajemng KN) ee catia) cay Se) gy Wat ke ecg rw mesiny knee Gre Te 


FROGRAM : MAXIMUM 


1 t 
! ' 
1 I 
' 
i \ 
t 
i THIS FROGRAM IS DESIGNED TO READ THE SQUARED ! 
: WAVEFORM VOLTAGE DATA AND TIME DATA STORED 
Lee St ove het CHOOSE. THE LOECSE Max Laur 
' (WAVEFORM ENSEMBLE? 2 THE CHUOSED MAXIMEMS WIEL 
: BE TAFEN SUARE-ROOT AND STORE THE FINAL RESULT , 
: Shine OLS 
t ' 
' 
' 


809 
810 
ea 
aa 
817 
814 
Slo 
216 j 
Bao DIM VDATAFlC20] ,VTIMESC20) .VNUMBERS£L 20] .VTIME MAXS£20) , VMAXIMUMSEL2O] 
900 DISF "WHAT IS THE NAME FOR VOLTAGE DATA FILE 7" 

910 INFUT VDATA 

20) DISF "WHAT IS THE NAME FOR TIME DATA FILE 7°" 

220 INFUT VTIMEE® 

940 DISF "WHAT IS THE NAME FOR FILE TO STORE THE NUMBER? NO. 7" 

950 INFUT VNUMBERS 

S60 PYSr BRAT 1S THE NAME FOR FILE TO STORE THE TIME MAX? DATA =" 
970 INFUT VTIME MAX 

980 DISF "WHAT IS THE NAME FOR FILE TO STORE THE MAXIMUMS?7" 

990 INFUT VMAXIMUNSS 

= DISF "WHAT IS THE STARTING TIME OF YOUR DATA COLLECING 7” 

992 INFUT START 

aay DISF "WHAT IS THE ENDING TIME OF THE DATA COLLECTING 7?" 

994 INFUT ENDT 

S96 

1900 ASSIGN# 1 TO VDATAS 

10929 ASSIGN# 3 TO VTIMEF 

1940 CREATE VNUMBER?,1,8 

19690 ASSIGN# 4 TO VNUMBERS 

ioe “SReArTE VIIME MAX? ,100,8 

1190 ASSIGN# 5 TO VTIME_MAK+ 

1120 CREATE VMAXIMNUNSS,100,8 

1140 ASSIGN# & TQ VMAXIMUMSF 

1141 : 

1144 INTEGER FRINTTER 

1146 DISF "WHAT IS THE FRINTER DEVICE ADDRESS NO. 7" 

1148 #! INFUT THE FPRINTER NO. SCREEN IS 1 

PLSoO INFUT FRINTTEF 

1152 FRINTER IS FRINTTER 

1154 : 

Vise Sores tisoly.1 (1501) ,TINC100) ,MAXS (100) 

ees Semin USING L170 + “TIME (SEC.>"., “VOLT” 

Lee IMAGE 25X%,11A,12X%,4A 

1180 INTEGER I,J 

1200 FOR I=START TO ENDT 

t229 READS 3 = TCI) 

1240 READ# 1 3: SCI) 

1260 NEXT I 

1280 LET J=9 

Peo 

b29 1 ! IF THE FIRST DATA BEING READ IS GREAT THAN THE SECOND ONE, 


Po en ee OOS ee See eee erate mec ise Sell fet Cen DON Ua Mien (oma Cen ) | © Egil. MO MEwiet cn) turache Ware Ureleaweencts teh Pe Mel ia, fei ie tay ey aie fae heh sia fe wo LWA BS '8) Paya Ree 
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a> > 
raed 
1.308 
Le) 
1340 
1360 
1780 
1400 
1410 
1420 
1440 
1460 
1480 
L500 
1520 
1S40 
eres) 
1580 
16400 
1510 
1620 
1440 
1460 
L680 
1700 
er ee 
f750 
1740 
1760 
1BoOo 
1829 
1840 
1880 
1900 


' THE FIRST OME IS TREATED AS 4 LOCAL Maximum. 
1 
IF S(START) “= S(START+1) THEN 1329 ELSE 1500 
TIM(G) =T START) 
MAXS(J)=SQR (S(START)) 
FRINT# S 3 TIM(I) 
FRINT# & 3: MAXS(J) 
FRINT USING 1410 3 "J=";9.TIM(J) .MAXS (J) 
IMAGE 12X,7A,DDD,7X,D.DDDDDD.11X.D. DDDDDDDD 
' PRINT TIM(J) 
| PRINT MAXS(J) 
J=J+1 
1 
FOR I=START+1 TO ENDT-1 
IF S(1).= $(I-1) AND S(I)2= S(I+1) THEN 1540 ELSE 1720 
VW Gait? 
MAXS(J)=SOR (S(1)) 
1 
FRINT USING 1610 3; “J=".Jd.TIM(J) MAXS(J) 
IMAGE 12X,54,DDD.7X,D.DDDDDD,11*%,D.DDDDPLDD 
! FRINT TIM(J) 
' FRINT MAXS(d) 
FRINT# S&S 3; TIM(J) 
FRINTH & : MAXS(J) 
J=eJg+ti 
NEXT I 
J=J-1 
FRINTH 4; J 
ASSIGNH 4 TO 
ASSIGN# 5 TO 
ASSIGN# & TO 
ASSIGN# 1 TO 
ASSIGN# 3 TO 
END 


Kok KK * 
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9oO0 
Gao} 
903 
904 
Sie 
906 
AZ. 
928 
Go9 
910 
<oled 
Sie 
oT 
914 
eles 
oe 
a1. 
mie 
oy 
1000 
Lo2o 
hie 
1924 
Leo 
10246 
1028 
19040 
1041 
19042 
1043 
1960 
LOBO 
1100 
eet) 
1140 
1146 
PIBoO 
1200 
1220 
1240 
$241 
1242 
1243 
1244 
12435 
L250 
Pest 
Zoe 


ea 
aU 


beaoe 
Za 
leo 
t2a% 
1260 
12860 
1200 
1320 
1340 
1360 


ee aa RIO er tae ese Alea ae I eae eS tet 


ee ee ee eee ee lees ee os Ce eae | OF me wi me) epee MNase” Geo inp im) wi Je) lay Ce) pe, "89 o snes ce Sco on eee ee 


EROGER Alec eOt. 


' | 
1 ' 
| 1 
| ' 
l 
: oe onhmh 1S DESt SHED OePeei hE ENV ELOE 
: OF WAVEFORIt BY READING DATA FROM DATA DIS. 
AFTER READING, VOLTAGE DATA IS CONVERTED TO LOG ; 
: PoOlCUE ter) SE PLOTTED, “AFTER FPLOMENG. A LEAST— 
} SGUBRE METHOR [5 USED TG FIND THE GEST FIT STRAIGHT . 
Pe ek sje Eee —-ENSENEEE PEO) THE LINE. ERINT 
; GUT THE SLOPE OF “THE JIVE Ane THE CORRELATION 
: SeQC-EITeClIeNT. : 
1 | 
J ! 
! J 


ee ee ee ee eae ee ee ee A ae I 


t 
t 
i 
| 
1 
t 
1 
' 
' 
t 
1 
I 
i 
1 
1 
ee ee oe ee er ee ee A tee ee Ae te 
: 

: 

CLEAR 

GCLEAF 

INTEGER FLOTER 

DISF “WHAT IS YOUF FOLTTEFE DEVICE ADDRESS NO. 3" 

‘ Me ouietHe FLOTTE Noe. SCREEM IS 1 

INFUT FLOTER 

FLOTTER IS FLOTER 

FRAME 

t 

THIS FLOT SIZE IS TO FIT THE REQUIREMENT FOR THE NFS THESIS 
t 

VOGETE 350,125,209 .,85 

CS i Z2e 

nove SO520 

ese Tine (“MICRO-SEL)”™ 

FEN UF 

MOVE 186,40 

DEG 

LDL 90 

ESBBEW “Lh ( VOLTAGE 2+) "P7so” 

FEN UF 

1 

| DEFAULT SET-UF IS TO FIX THE FLOTING COORDINATION AS 
X-AXIS IS FROM Son TO 670 (MICRO-SEC), Y-AXIS IS FROM 

! -§ TO 2 

' 

Dise DO yOu WANT FO USE DEFAULT SET-UP = YES=0" 

INFUT BULL 

Deaesuek=o THEN 1255 ELSE 1260 

Sea U0, 7 so. = 

Seen Moo Jeon sae 

AMIN=S500 

GOTO 1710 

1 

DiS. ENTER THe Xidiil! OF SCALE” 

INFUT XMIM 

DISF "ENTER XMAX OF SCALE” 

INFUT XMAX 

Piste eigerm THE YMIN OF SCALE” 

INFUT YMIN 


ql 


. 


1280 
14009 
1420 
1440 
1460 
1489 
1500 
1320 
1540 
15460 
1580 
1400 
146420 
1649 
1669 
14661 
1480 
1709 
Lye 
L772 
17 27 
1740 
1741 
1742 
1744 
1746 
1748 
1749 
i7sao 
eo 
Ie 
lef 
17460 
1780 
17H 
176e 
1730) 
1800 
1802 
1803 
1804 
1810 
1820 
1821 
leee 
1e235 
Legs 
hee 
1828 
te. 


LBS 


1823 
1834 
1835 
18737 
1839 
1840 
1726) 
1930 
1940 


DISF "ENTER THE YMAX CF SEALE ® 

INFUT YMA 

SCALE XMIN,XMAX,YMIN,YMAX 

XD ee 

DISF “ENTER THE X=-TICr ING SFACE" 

INPUT 

DISF "ENTER THE Y-TICt ING SFACE" 

INPUTS yt 

DISF "ENTER THE X INTERSECTION" 

INP ERAa 

DISF EM Ger rey eee oe lr. 

INPUTS ya 

DISF "ENTER THE X=MAJOR COUNT" 

INFUT XME 

DISF "ENTER THE Y-MAJOR COUNT" 

' 

INFUT YHCe 

LAXES XU yee leet leer, nies 

GRAFH 

Bre iNael 

1 

DIM- TIN CiIOO}. MARS (190? 

DIM” TInEMAK2£ 200). Mak TMS beng 

DISF "WHAT IS YOUR FILE NAME FOF TIMEMAX? 7" 
INFUT TIMEMAX = 

DISF "WHAT IS YOUR FILE NAME FOR MAXIMUM: 7" 
INPUT MAX IMUMS 

FRINT "THE DATA FILES ARE" ,TIMEMAX? , "AND" ,MAxX IMUMs 
DISF “WHAT IS THE HIGHEST ORDER OF AREAY NUMER” 
' THIS GIVES A CHOOICE TO CHOOSE THE END FOINT OF FLOT 
INFUT JK 

| 

ASSIGN# 5 TO TIMEMAX+S 

ASSIGN# 6&6 TO MAXIMUMS 

DISF “FRESS © IF THE FIRST LOCAL MAXIMUM IS BEING SAVED 7 
INPUT OF 

FOR I=0 TO JE 

READ# S&S : TIM(CI) 

{ 

‘ CONVERT THE TIME UNIT FROM SEC le shileRags == 
l 

TIMCIV=TIM( I) *¥1000000 

READH & =: MAXS(I) 

MAXS(I)=LOG (MAXS(I)) 

: CONVERT THE VOLTAGE TO NATURE LOG VALUE 

NEXT I 

K=Jt+1 

IF OKk=0 THEN GOTO i835 

FOR a=1 Todi 

The 1) = eric) 

MAXS(I-1)=MAXS(I) 

NEXT I 

K=JE. 

' DISF “WHICH CHARACTER YOU FREFER TO USING IN FLOTTING THE 
' INFUT C# 

MOVE TIM(O) ,MAXS(0) 

FOR I=0 TO k-1 

DRAW TIMCI) .MAXS(I1) 

' LABEL CF 

NEXT I 
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DATAe: 


rou 
LS | 
to. 
1980 
2QO00 
at! 
2900 
=7 10) 
292C 
ae 
29s4 
So 
=726 
27a) 
SOOO 
2040 
3080 
TOBO 
3100 
5120 
3140 
alow 
S140 
SLO 
Soo 
S240 
S20 
2280 
ooo 
Sane) 
30340 


pe ee) 


BS 


Soe) 
SS 
Soe) 
S570 
Soe 1 
Soy 
3580 
3400 
32420 
3421 
547f 
3422 
35440 
2441 
3442 
2443 
=460 
3461 
3480 
3481 
SSO0 
Saal 
ce 
2 2wI4AOD 
3560 
5a) | 
ao 


voaeO 


er es 

t 

GOSUS LEASTSCR 
ASSIGN# 3S TO * 
ASSIGN# 6 TO * 
END 


FGELGWING SUBROUTINE [S DESIGNED TO APFLY THE LEAST 
SQUARE METHOD IN THE ABOVE DATA AND THE FLOT THE 
BEST Fit StRaIGHT EINE. FRINT OUT THE SLOFE AND 

(ies eCOKRelat [OM COSEEIEIENT ON THE FLOT. 

Ponnlirm USED FEEASE 


n 
rm 


Seek: 
AL Y,A4Y,X%1,X£,TAVE,YAVE,VARX , VARY ,COVXY 
T Y=0 

Y I5 THE SUMMATION OF YCI) VALUE 
XY=0 

XY IS THE SUMMATIGN OF X(T) *¥CT1) 
KX 1=0 

Meloy idee olnmriiimni. ON OF x Cl) 


e 
1S 
cs 


aii a (= a ees ee oe = 


- =o 
in om 
4 


EST. 2-6 

: Fool Sate. SUMMA) ONT OF 2s CL)ooS 
ROR ONG k-T1 

Y=Y+MAXS (1) 

XY=AY+TIMCI) «MAXS CI) 

X1=X1+TIM(1) 

Re=AS ei itt 1) 2 

NE x ined 

B= CX ex —-% PHXY) / (Ch eX S—X PQ) 


: AopltometHe SLOPE OF THE PItlee SikAlGH Se rent 
Gara KY-ALAY) / (Rh *#XS-X1° 2) 
i 
' 
i 


Beets THe INTEREEPT TO4y 3440S 
Poltiestne OLOPE A=" sr 

Seine y—INTEREEPT Baars 
YAVE=(MAXS (0) +MAXS (R-1)) 72 


YAVE IS THE AVERAGE VALUE FOR THE LOG VOLTAGE 
AVE=(TIM(O0)4+TIM(EK-1))7/2 


TAVE IS THE AVERAGE VALUE FOR THE TIME VARIABLE 


COVXY IS THE COVARIANCE OF TIME AND VOLTAGE (LN) 


VARX=0O 

' VARX IS THE VARIANCE OF TIME VARIABLE 

b 

VARY=0 

t VARY IS THE VARIANCE OF VOLTAGE(LN) VARIABLE 
4 

DIM YEC100) 

FOR I=0 TO k-1 


ds 


£400 
S620 
A400 
25640 
24689 
e700 
Jeo 
3740 
3740 
3741 
2780 
2B00 
Be = 
3840 
2860 
7~BB8O 
Z9Q0 
ee) 
S720 
2740 
ao) 
39789 
eae 
299 1 
a924 
Seo 
Zoo 
O79 9 
4oog 
4002 
4OD4 
4008 
4020 
4921 
4n40 
4060 
4061 
4042 
4080 


5 

COVXY=COVXY+(TIM( 1) -TAVE) * (MAXS (1) -YAVE) 

5 

VAR X=VARX+ (TIM(1) -TAVE) 72 

t 

VARY=VARY + (MAXS (1) -YAVE) °2 

1 

YL(1)=A*TIM(1) +B 

! YL IS THE VALUE FITTED ON THE STRAIGHT LINE 
| 

NEXT I 

FEN 2 

MOVE TIM(O) ,YL (9) 

FOR I=0 TO K-1 

DRAW TIM(T) ,YL(T) 

NEXT I 

CORCOE=COVXY/SOQR (VARX#*VARY) 

t 

! CORCOE IS THE CORRELATION COEFFICIENT 

4 

PRINT “CORRELATION COEFFICIENT IS ";CORCOE 

FEN 1 

DEG 

LDIR 9 

4 

FOLLOWING RRINT OUT THE SLOFE &ND COORELATION COEFFICIENT 
ON THE FLOT 

| 

MOVE XMIN+80, 1 

' 

LABLE ALWAYS STARTS FROM THE XMIN+89 IN X COORDINATES 
| 

LABEL USING 4021 ; "SLOFE =",A 

IMAGE 7A,DDD.DDDD 

MOVE XMIN+80,2 

LABEL USING 40461 ; “CORRELATION COEFF=",CORCOE 
IMAGE 19A,DDD.DDDD 

FEN UF 

RETURN 


74 


Boy 
Bi 
B20 
Gz 
840 
ye 
860 
870 
B80 
8°90 
1900 
1020 
1046 
1060 
1980 
1109 
1120 
1140 
11460 
1180 
1200 
P2209) 
Pace 
ee 


Sa 


232 
1240 
tae) 
1261 
L2G 
1263 
127 @ 
27 1 
Tere 
1274 
ie 
1276 
1230 
13500 
cet 
1202 
1340 
126O 
Lay 1 
1S72 
PSs 
1580 
1400 
142 
1440 
Lo00 
1520 
1540 


THE DATA SETS TO GET A AVERAGED 


: 

! 

FROGRAM : AVERAGE 
1 

' 

| 

: WAVEFORM DATA. 


1 
J 
' 
' 
THIS PROGRAM 1S DESIGNED TO AVERAGE 
I 
| 
' 
{ 


. . . 


eee eee Oe net ee. ee a6 e716) el 6 (est ete 6) ie #) a 6. @ (i "e je «el (ie. (se eo, 6) oo 3 le be je 16 60 er *et o> Se (@) Oe 8: 


DIM Vi (1500) ,Vv2(1500) ,v3(1500) 

DISF "WHAT IS THE ORIGINAL DATA FILE 7? DAT#H" 
INFUT VDAT# 

DISP “WHAT IS THE NEW DATA FILE 7 DATA#" 
INPUT VDATAS 

DISP “WHAT IS THE AVERAGED DATA FILE 7 DATi#" 
INFUT VDATI4 

ASSIGN# 1 TO VDATS 

ASSIGNH 2 TO VDATAF 

CREATE VDAT1#,8,110 

ASSIGN# 3 TO vDaTiF 

i 

! N IS THE NUMBER HOW MANY DATA SET BEING AVERGED 
| 

DISF “WHAT IS YOUR AVERAGE No. N 7?" 

INFUT N 


' AVERAGE ONLY AFFLY TO THE SAME ARRAY NO. 


mul l=s22 10 621 
1 


READ# pemew) 4 1) 
NEA et 

PORew a2? 10 621 
READH 2 3; Val) 
NEXT I 


POR l=a22 10 621 
RUNNING AVERAGE IS ADOFTED 


Var NF) YI CI Ve CT) ZN 
PFRINT# 3 3: VSCt) 


' PRINT OUT FOLLOWING VALUE AS A CHECE ING FURFOSE 
1 

IF T=S522 THEN 1500 ELSE 1400 

IF JT=55S THEN 1500 ELSE 1420 

IF I=600 THEN 1500 ELSE 1440 

IF I=621 THEN 1500 ELSE 1520 

PRibiiaviti >) Vv2tl) ,Vs(l) 

NEXT I 

END 


V2 


APPENDAXK C 
DECAY ENVELOPE GRAPHS FOR AGGREGATE 


This appendix includes the graphs of the decay envelopes for reflection from 
aggregate water interface within the glass tank. Figure C.1 to C.9 show the original 
decay envelopes obtained from different transducer positions. The decav envelopes are 
also described as following data sets: 

data set 42 
data set 43 
data set 44 
data set 45 
data set 46 
data set 47 
data set 48 
data set 49 
data set 50 
Figures C.10 to C.18 are the decay envelopes resulted from the running average 
method. They are also described as following data sets: 
data set DAT 2 
data set DAT3 
data set DAT4 
data set DAT5 
data set DAT6 
data set DAT7 
data set DATS 
data set DAT9 
data set DAT10 
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LN ((<VOLTAGE*2>) 71/2) 
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Figure C.1 Decay envelope for reflection from aggregate/water interface, 
(ralisducen eas SU cI above Slillace (data set 42). 
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figure C.2 Decay envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 43). 
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Figure C.3 Decay envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 44). 
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Pigure C.4 Decay envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 45). 
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Figure C.5 Decay envelope for reflection from a 
transducer was 30 cm above surface ( 
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Figure C.6 Decay envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 47). 
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Figure C.7 Decav envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 48). 
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Figure C.8 Decay envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 49). 
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Pigure C.9 


Decav envelope for reflection from aggregate/water interface, 
transducer was 30 cm above surface (data set 5 
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envelope resulted from the average of data set 42 and 43 
agercgate/water surface (data set DAT 2). 
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Figure C.11 
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Decay envelope resulted from the average from data set 42 to 44 
for the aggregate/water surface (data set DAT 3). 
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Figure C.12) Decay envelope resulted from the average from data set42 to 45 
for the aggregate/water surface (data set DAT4). 
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Figure C.13) Decay envelope resulted from the average from data set 42 to 46 
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for the aggregate/water surface (data set DAT5). 
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Figure C.14 easy envelope resulted from the average yom data Set Ao 


for the aggregate/water surface (data set DATO). 
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Figure C.15 Decay envelope resulted from the average from data set 42 to 48 
for the ageregate/water surface (data set DAI 7). 
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Figure C.16 Decav envelope resulted from the average from data set 4? to 49 
or the aggregate/water surface (data set DATS). 
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Figure C.17) Decay envelope resulted from the average from data set 42 to 50 
for the aggregate/water surface (data set DA 19). 
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Figure C.18 Decay envelope resulted from the av ae Ecotn jy set 42 to 5! 
for the aggregate/water surface (data set DATI 
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APPENDIX D 
Weer PiNVELOPE GRAPHS FOR FINE SAND 


This appendix includes the graphs of the decay envelopes for reflection from fine- 
sand, water interface within the glass tank. Figure D.1 to D.10 show the original decay 
envelopes obtained from different transducer positions. The decay envelopes are also 
described as following data sets: 

Gala SCL Oo 
data set 56 
data set 57 
data set 58 
data set 59 
data set 60 
data set 6! 
data set 62 
data set 63 
data set 64 

Figures D.11 to D.19 are the decay envelopes resulted from the running average 

method. They are also described as following data sets: 
data set DAT12 
data set DAII3 
data set DATI4 
data set DATI5 
data ser WAT 1G 
data set DATI7 
data set DATI8 
data set DATI9 
data set DAT20 
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Figure D.1 Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 55). 
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Figure D.2. Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 56). 
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Decav envelope for reflection from fine-sand/water interface, 


transducer was 30 cm above surface (data set 57). 
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Figure D.4 Decay envelope for reflection from fine-sand; water interface, 
transducer was 30 cm above surface (data set 58). 


oe 


Lit Gia Ge ee) 17 2) 





f 
oe 





= oe 
[T] ’ 
ee | 
a f 
o { 
© 560). 
in | 
T 
O 

Sour ep Go 

Gr 0 

: 2 OB 

| ’ 

\ l sie 

: = 

GOS f¢- 1 5 

, a 

r = c 

a Oo m7 

! T7 

620 > A 

} 1 

2 

640 f- S) 
r 
BGC f+ 


Figure D.5 Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 59). 
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Figure 10.6 Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 6Q). 
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Figure D.7 Decay envelope for reflection from fine-sand/water interface, 
transducer was 30 cm above surface (data set 61). 
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Decay envelope for reflection from fine-sand‘ water interface, 


transducer was 30 cm above surface (data set 62). 
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Figure D.9 Decav envelope for reflection from fine-sand,water interface, 
transducer was 30 cm above surface (data set 63). 
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Figure D.10 Decay envelope for reflection from fine-sand/ water interface, 
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transducer was 30 cm above surface (data set 64). 
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Figure D.11 Decay envelope resulted from the average of data set 55 and 56 
for the fine-sand/water surface (data set DAT12). 
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Figure D.12 IWecay envelope resulted from the average from_data set 55 to 57 
or the fine-sand/water surface (data set DAT 13). 
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Figure D.13 > Decay envelope resulted from the average from data set 55 to 58 
or the fine-sand/water surface (data set DAT 14). 
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Figure D.14 Decay envelope resulted from the average from data set 55 to 59. 
or the fine-sand/water surface (data set DATI5). 
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Figure D.15 Decay envelope resulted from the average from data set 55 to 60 
or the fine-sand/water surface (data set DAT 16). 
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Figure D.16 Decay envelope resulted from the average from_data set 55 to 61 
for the f{ine-sand/water surface (data set DAT 17). 
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Figure D.17 Decay envelope resulted fropitliior ena casino imc atdesctsopmtiarcy 
for the fine-sand/water surface (data set DAT 18). 
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Decay envelope resulted from the average from data set 55 to 63 
for the fine-sand, water surface (data set DATIY). 





LNG (C<VOLTAGE* 2>) 71/2) 


| { 1 
Ci + G) rv 


920 


ie 


940 


= 


eon) 


& 560 


c 
Od end end 


i] 


080 


600 


640}- 


660 


Pigure D.19 Decay envelope resulted from the average from data set 55 to 64 
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for the aggregate/water surface (data set DAT 20). 
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